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We calculate the impact of the complete set of two-pion exchange contributions at chiral order 
four (also known as next-to-next-to-next-to-leading order, N 3 LO) on peripheral partial waves of 
nucleon-nucleon scattering. Our calculations are based upon the analytical studies by Kaiser. It 
turns out that the contribution of order four is substantially smaller than the one of order three, 
indicating convergence of the chiral expansion. We compare the prediction from chiral pion-exchange 
with the corresponding one from conventional meson-theory as represented by the Bonn Full Model 
and find, in general, good agreement. Our calculations provide a sound basis for investigating the 
issue whether the low-energy constants determined from irN lead to reasonable predictions for NN. 
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O ■ I. INTRODUCTION 
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■ One of the most fundamental problems of nuclear physics is to derive the force between two nucleons from first 
rS^ ! principles. A great obstacle for the solution of this problem has been the fact that the fundamental theory of strong 
'"^ ■ interaction, QCD, is nonperturbative in the low-energy regime characteristic for nuclear physics. The way out of this 

dilemma is paved by the effective field theory concept which recognizes different energy scales in nature. Below the 
chiral symmetry breaking scale, A x « 1 GeV, the appropriate degrees of freedom are pions and nucleons interacting 
via a force that is governed by the symmetries of QCD, particularly, (broken) chiral symmetry. 

The derivation of the nuclear force from chiral effective field theory was initiated by Weinberg and pioneered by 
Ordonez Q and van Kolck Subsequently, many researchers became interested in the field ||5|-|20|] . As a result, 

(*C) ', efficient methods for deriving the nuclear force from chiral Lagrangians emerged 0H^] and the quantitative nature 

■ of the chiral nucleon-nucleon (NN) potential improved fll^O]. 

\ Current JVJV potentials Jl^JlJl and phase shift analyses [)21|| include 27r-exchange contributions up to order three 
in small momenta (next-to- next-to-leading order, NNLO). However, the contribution at order three is very large, 
several times the one at order two (NLO). This fact raises serious questions concerning the convergence of the chiral 



> 



expansion for the two-nucleon problem. Moreover, it was shown in Ref. |14| that a quantitative chiral NN potential 
' requires contact terms of order four. Consistency then implies that also 2tt (and 3n) contributions are to be included 
I , up to order four. 

For the reasons discussed, it is a timely project to investigate the chiral 2tt exchange contribution to the JVJV 
interaction at order four. Recently, Kaiser pd| , ^2| has derived the analytic expressions at this order using covariant 
■ perturbation theory and dimensional regularization. It is the chief purpose of this paper to apply these contributions 
in peripheral NN scattering and compare the predictions to empirical phase shifts as well as to the results from 
conventional meson theory. Furthermore, we will investigate the above-mentioned convergence issue. Our calculations 
provide a sound basis to discuss the question whether the low-energy constants (LECs) determined from 7rJV lead to 
reasonable predictions in NN. 

In Sec. II, we summarize the Lagrangians involved in the evaluation of the 27r-exchange contributions presented 
in Sec. III. In Sec. IV, we explain how we calculate the phase shifts for peripheral partial waves and present results. 
Sec. V concludes the paper. 

II. EFFECTIVE CHIRAL LAGRANGIANS 

The effective chiral Lagrangian relevant to our problem can be written as |22 23 1 , 

r _ r (2) , /.(l) , r (2) , r {3) , /-.x 
*-cff — inr T *-7rjV ' '"irN 7T N T • • • , 
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where the superscript refers to the number of derivatives or pion mass insertions (chiral dimension) and the ellipsis 
stands for terms of chiral order four or higher. 

At lowest/leading order, the irir Lagrangian is given by, 



4 2 i = f -f tr [d»Ud^ + ml(U + t/t) ] 



(2) 



and the relativistic irN Lagrangian reads, 
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The coefficient a that appears in the last equation is arbitrary. Therefore, diagrams with chiral vertices that involve 
three or four pions must always be grouped together such that the a-dependence drops out (cf. Fig. |[ below). 

In the above equations, Mn denotes the nucleon mass, gA the axial- vector coupling constant, and f n the pion decay 
constant. Numerical values are given in Table I. 

We apply the heavy baryon (HB) formulation of chiral perturbation theory p3] in which the relativistic ttN 
Lagrangian is subjected to an expansion in terms of powers of 1/Mjy (kind of a nonrelativistic expansion), the lowest 
order of which is 
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In the relativistic formulation, the field operators representing nucleons, 'J, contain four-component Dirac spinors; 
while in the HB version, the field operators, N, contain Pauli spinors; in addition, all nucleon field operators contain 
Pauli spinors describing the isospin of the nucleon. 

At dimension two, the relativistic irN Lagrangian reads 



4^ = E^ 2) *- 



(9) 



The various operators Oj 2 are given in Ref. p3| . The fundamental rule by which this Lagrangian — as well as all the 
other ones — are assembled is that they must contain all terms consistent with chiral symmetry and Lorentz invariance 
(apart from the other trivial symmetries) at a given chiral dimension (here: order two). The parameters c; are known 
as low-enery constants (LECs) and are determined empirically from fits to ttN data (Table I). 
The HB projected irN Lagrangian at order two is most conveniently broken up into two pieces, 



5(2) _ £(2) 
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Note that C^L fi x is created entirely from the HB expansion of the relativistic C^jj and thus has no free parameters 
2) 

("fixed"), while C^jq ct is dominated by the new nN contact terms proportional to the parameters, besides some 
small 1/Mtv corrections. 

At dimension three, the relativistic irN Lagrangian can be formally written as 

23 
i=l 

with the operators, 0\ , listed in Refs. p2| , p3[ ; not all 23 terms are of interest here. The new LECs that occur at 
this order are the di. Similar to the order two case, the HB projected Lagrangian at order three can be broken into 
two pieces, 

ri 3 ) _ ri 3 ) I r( 3 ) (-ia) 

7riV — '-ttN, fix ' ^TrJV.ct' V 1 ^ 



with ^Jfjv fi x and £^ ct given in Refs. p2p3|. 



III. NONITERATIVE 2tt EXCHANGE CONTRIBUTIONS TO THE NN INTERACTION 

The effective Lagrangian presented in the previous section is the crucial ingredient for the evaluation of the pion- 
exchange contributions to the nucleon-nucleon (NN) interaction. Since we are dealing here with a low-energy effective 
theory, it is appropriate to analyze the diagrams in terms of powers of small momenta: (Q/A x ) 1 ' , where Q stands 
for a momentum (nucleon three-momentum or pion four-momentum) or a pion mass and A x ps 1 GeV is the chiral 
symmetry breaking scale. This procedure has become known as power counting. For non-iterative contributions to 
the NN interaction (i. e., irreducible graphs with four external nucleon legs), the power v of a diagram is given by 

^ = 21+^(^ + ^-2) , (15) 
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where I denotes the number of loops in the diagram, dj the number of derivatives or pion-mass insertions and 
the number of nucleon fields involved in vertex j; the sum runs over all vertices j contained in the diagram under 
consideration. 

At order zero (u = 0, lowest order, leading order, LO), we have only the static one-pion-exchange (OPE) and, at 
order one, there are no pion-exchange contributions. Higher order graphs are shown in Figs. |l|-|^. Analytic results for 
these graphs were derived by Kaiser and coworkers [[7 pTl]Jl2] | using covariant perturbation, i. e., they start out with 
the relativistic versions of the irN Lagrangians (see previous section). Relativistic vertices and nucleon propagators 
are then expanded in powers of 1/Mjv- The divergences that occur in conjunction with the four-dimensional loop 
integrals are treated by means of dimensional regularization, a prescription which is consistent with chiral symmetry 
and power counting. The results derived in this way are the same obtained when starting right away with the HB 
versions of the irN Lagrangians. However, as it turns out, the method used by the Munich group is more efficient in 
dealing with the rather tedious calculations. 

We will state the analytical results in terms of contributions to the on-shell momentum-space NN amplitude which 
has the general form, 

V(p',p) = { V c +t 1 -t 2 W c 



(27T) 



+ [Vs + T% ■ T 2 Ws ] C?l • <?2 

+ [ Vt + Ti ■ r 2 Wt ] <?i • q <?2 • q 
+ [V ls + t 1 - t 2 W LS ] i-iS -(qxk) 



+ [V aL +r 1 ■T 2 W aL ] ai-{qxk) a 2 -{qxk)j, (16) 

where p ' and p denote the final and initial nucleon momentum in the center-of-mass (CM) frame, respectively, 

q = p ' — p is the momentum transfer, 
k = i (p ' + p) the average momentum, 
£ = ^((?i + a 2 ) the total spin, 
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and CT1.2 and n 2 are the spin and isospin operators, respectively, of nucleon 1 and 2. For on-energy-shell scattering, 
V a and W a (a — C, S, T, LS", trL) can be expressed as functions of q and k (with q = \q\ and fc = |fc|), only. 

Our formalism is similar to the one used in Refs. |?|-[l2]|, except for two differences: all our momentum space 
amplitudes differ by an over-all factor of (—1) and our spin-orbit amplitudes, Vls ancl Wls, differ by an additional 
factor of (—2) from the conventions used by Kaiser et al. 
closely in tune with what is commonly used in nuclear physics. 

We stress that, throughout this paper, we consider on-shell NN amplitudes, i. e., we always assume \p'\ = \p\ = p. 
Note also that we will state only the nonpolynomial part of the amplitudes. Polynomial terms can be absorbed into 
contact interactions that are not the subject of this study. Moreover, in Sec. IV, below, we will show results for NN 
scattring in F and higher partial waves (orbital angular momentum L > 3) where polynomials of order Q v with v < 4 
do not contribute. 



- 1 2 1 . We have chosen our conventions such that they are 



A. Order two 



Two-pion exchange occurs first at order two [y = 2, next-to-leading order, NLO), also know as leading-order 2ir 
exchange. The graphs are shown in the first row of Fig. [|. Since a loop creates already v = 2, the vertices involved 
at this order can only be from the leading/lowest order Lagrangian C^ Nl Eq. (||), i. e., they carry only one derivative. 
These vertices are denoted by small dots in Fig. [j]. Note that, here, we include only the non- iterative part of the box 
diagram which is obtained by subtracting the iterated OPE contribution [Eq. J6q), below, but using Mf^/E p w Mjv] 
from the full box diagram at order two. To make this paper selfcontained and to uniquely define the contributions 
for which we will show results in Sec. IV, below, we summarize here the explicit mathematical expressions derived in 
Rcf. @: 



W c = 
V T = 



384vr 2 /^ 



~ Vs = 
q 



4ml(5g 4 A - Ag\ - 1) + q 2 {2ig\ - U)g\ - 1) 
64tt 2 ^ ' 



48g%m* 



(17) 
(18) 



where 



L{q) 



w w + q 
— m 

q Im-p; 



(19) 



and 



\J 4m 2 + q 2 



(20) 



B. Order three 



The two-pion exchange diagrams of order three (y = 3, next-to-next-to-leading order, NNLO) are very similar to the 
ones of order two, except that they contain one insertion from Eq. (|l0|). The resulting contributions are typically 
either proportional to one of the low-energy constants Cj or they contain a factor 1/Mjy. Notice that relativistic 1/Mjv 
corrections can occur for vertices and nucleon propagators. In Fig. ^, we show in row two the diagrams with vertices 
proportional to (large solid dot), Eq. Jl2]), and in row three and four a few representative graphs with a 1/M/v 
correction (symbols with an open circle). The number of 1/Mjv correction graphs is large and not all are shown in 
the figure. Again, the box diagram is corrected for a contribution from the iterated OPE: in Eq. j65|), below, the 
expansion of the factor Al^/E p — Mjv — P 2 /Mn + ■ ■ ■ is applied; the term proportional to (—p 2 /Mjy) is subtracted 
from the third order box diagram contribution. For completeness, we recall here the mathematical expressions derived 
in Ref. §. 



V c = 
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9\ 
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with 
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C. Order four 



This order, which may also be denoted by next-to-next-to- next-to-leading order (N 3 LO), is the main focus of this 
paper. There are one- loop graphs (Fig. |J) and two- loop contributions (Fig. |3|). 



1. One-loop diagrams 



a. c 2 contributions The only contribution of this kind comes from the football diagram with both vertices pro- 



portional to Ci (first row of Fig. |2j). One obtains |ll| 

3L(q) 



V c 



W T = — ~W S 



—w + c 3 w - 4cim w J + —w 
c\w 2 L(q) 

967T2/4 • 



(29) 
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b. Ci /Mn contributions This class consists of diagrams with one vertex proportional to c; and one 1 /-Mjy correction. 
A few graphs that are representative for this class are shown in the second row of Fig. |[ Symbols with a large solid 
dot and an open circle denote 1/Mjv corrections of vertices proportional to c,. They are part of cfh' ^q. (|l4|). The 
result for this group of diagrams is pi| : 



V c 
W c 
W T 
Vls 
W LS 



9 2 A L{g) 
32ir 2 M N f* 
CAq 2 L{q) 
' 192tt 2 M N f* 
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(31) 
(32) 
(33) 
(34) 
(35) 
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c. 1/Mjy corrections These are relativistic 1/Mjy corrections of the leading order 2ir exchange diagrams. Typical 
examples for this large class are shown in row three to six of Fig. ^. This time, there is no correction from the iterated 

{jt/Ep does not create a term proportional to 1/Mjy. 



OPE, Eq. (|65]), since the expansion of the factor M^/E p does not create a term proportional to 1/Mjy. The total 



result for this class is |L2l 
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In the above expressions, we have replaced the p 2 dependence used in Ref. |l2| by a k 2 dependence applying the 
(on-shell) identity 



2 1 2 i r 2 
p = -q +k . 



(43) 



2. Two-loop contributions 



The two-loop contributions are quite involved. In Fig. 0, we attempt a graphical representation of this class. The 
gray disk stands for all one-loop TtN graphs which are shown in some detail in the lower part of the figure. Not all of 
the numerous graphs are displayed. Some of the missing ones are obtained by permutation of the vertices along the 
nucleon line, others by inverting initial and final states. Vertices denoted by a small dot are from the leading order nN 
Lagrangian L^ N , Eq. (||), except for the 4ir vertices which are from C^l, Eq. (||). The solid square represents vertices 



proportional to the LECs di which are introduced by the third order Lagrangian Z^y, Eq. (13). The di vertices occur 
actually in one-loop NN diagrams, but we list them among the two-loop NN contributions because they are needed to 
absorb divergences generated by one-loop ttN graphs. Using techniques from dispersion theory, Kaiser jnj calculated 
the imaginary parts of the NN amplitudes, Im V a {in) and Im W a {i^), which result from analytic continuation to 
time-like momentum transfer q = ifi — + with /i > 2m n . We will first state these expressions and, then, further 
elaborate on them. 
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where k = ^//U 2 /4 — m 2 . 

We need the momentum space amplitudes V a (<7) and W a (<z) which can be obtained from the above expressions by 
means of the dispersion integrals: 
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and similarly for Wc,s,t- 

We have evaluated these dispersion integrals and obtain: 
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with 
and 



w } (9) = -^{q) = - y t;:^ (du d 15 ) (so 



^ = ~W S = — L T^WVq^ ' (61) 



Wr(q) = ~Ws(q) = ^0JW [™ 2 MV) + 2m 7r (l + 2<&)] (62) 

We were able to find analytic solutions for all dispersion integrals except W^f 1 and (and V^). The analytic 
solutions hold modulo polynomials. We have checked the importance of those contributions where the integrations 
have to be performed numerically. It turns out that the combined effect on NN phase shifts from W^f 1 , , and Vg b ^ 
is smaller than 0.1 deg in F and G waves and smaller than 0.01 deg in H waves, at 71 a b = 300 MeV (and less at lower 
energies). This renders these contributions negligible, a fact that may be of interest in future chiral NN potential 
developments where computing time could be an issue. We stress, however, that in all phase shift calculations of this 
paper (presented in Sec. IV, below) the contributions from Wq\ \ and are always included in all fourth order 
results. 

In Eqs. (|57|) and (|60|), we use the scale-independent LECs, d\, which are obtained by combining the scale-dependent 
ones, <^(A), with the chiral logarithmus, ln(m T /A), or equivalently dj = d^im^). The scale-dependent LECs, <i£(A), 
are a consequence of renormalization. For more details about this issue, see Ref. p2[ . 



IV. NN SCATTERING IN PERIPHERAL PARTIAL WAVES 



In this section, we will calculate the phase shifts that result from the NN amplitudes presented in the previous 
section and compare them to the empirical phase shifts as well as to the predictions from conventional meson theory. 
For this comparison to be realistic, we must also include the one-pion-exchange (OPE) amplitude and the iterated 
one-pion-exchange, which we will explain first. We then describe in detail how the phase shifts are calculated. Finally, 
we show phase parameters for F and higher partial waves and energies below 300 MeV. 



A. OPE and iterated OPE 



Throughout this paper, we consider neutron-proton (rip) scattering and take the charge-dependence of OPE due to 
pion-mass splitting into account, since it is appreciable. Introducing the definition, 

V w (m n ) = -Tp^-^ —2— - — — , 63) 

the charge-dependent OPE for np scattering is given by, 

V OPE (p',P) = -K(m 7r D) + (-l) /+1 2F 7r (m 7r ±), (64) 

where I denotes the isospin of the two-nucleon system. We use m,o = 134.9766 MeV and m 7r ± — 139.5702 MeV p6[ . 
The twice iterated OPE generates the iterative part of the 27r-exchange, which is 



v 2 „, it (p',p) = ^- J dV 



1 1 „ V ope (p',p")Vopk(p",P) 



p2 _ p/l _|_ j £ 

where, for Afjv, we use twice the reduced mass of proton and neutron, 

2M p M n 



(65) 



M N = p " = 938.9182 MeV , (66) 



M p + M, 



and E p ee ^fM 2 N +p 2 . 

The T-matrix considered in this study is, 
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T(p',p) = VoPE(p",p) + V 27r at(p',P) + V 27 r,irr(p',P) 



(67) 



where Vi-k^tt refers to any or all of the contributions presented in Sec. III. In the calculation of the latter contributions, 
we use the average pion mass = 138.039 MeV and, thus, neglect the charge-dependence due to pion-mass splitting. 
The charge-dependence that emerges from irreducible 2-k exchange was investigated in Ref. [^9| and found to be 
negligible for partial waves with L > 3. 



B. Calculating phase shifts 



We perform a partial-wave decomposition of the amplitude using the formalism of Refs. [|30|-|32[ . For this purpose, 
we first represent T(p',p), Eq. (|67|), in terms of helicity states yielding (p 'A^ \' 2 \T\fi\i\2). Note that the helicity A; 



of particle i (with i = 1 or 2) is the eigenvalue of the helicity operator \tii ■ Pi/\pi 
angular momentum states is accomplished by 

(X' 1 X' 2 \T J (p',p)\X 1 X 2 ) = 2n 



which is ±|. 



Decomposition into 



(68) 



where 9 is the angle between p ' and p and d J mm , (9) are the conventional reduced rotation matrices which can 
be expressed in terms of Legendre polynominals Pj (cos 9). Time-reversal invariance, parity conservation, and spin 
conservation (which is a consequence of isospin conservation and the Pauli principle) imply that only five of the 16 
helicity amplitudes are independent. For the five amplitudes, we choose the following set: 



T{(p,p) = (+ 


+ \T J (p,p)\ + +) 


T/(p,p) = (+ 


+ \T J (p,p)\--) 


T 3 J (p,p) = (+ 


-\T J (p,p)\ + -) 


T({p, P ) = (+ 


-\T J (p,p)\-+) 


T 5 J (p,p) = (+ 


+ \T J (p,p)\ + -) 



(69) 

where ± stands for ±-|, and where the repeated argument (p,p) stresses the fact that our consideration is restricted 
to the on-shell amplitude. The following linear combinations of helicity amplitudes will prove to be useful: 

°T J = T( - T$ 
i t j = t j _ t j 

12 T J = T J + T J 



3irpJ _ rpj 



2T( 



(70) 



More common in nuclear physics is the representation of two-nuclcon states in terms of an \LSJM) basis, where S 
denotes the total spin, L the total orbital angular momentum, and J the total angular momentum with projection 
M. In this basis, we will denote the T matrix elements by T/, 5 L = (L 1 SJM\T\LSJM). These are obtained from the 
helicity state matrix elements by the following unitary transformation: 
Spin singlet 



— 



Uncoupled spin triplet 
Coupled triplet states 



T 



.71 



J.J 



OrpJ 



IrpJ 



(71) 
(72) 



T 



Ji 

j-i.j-i 



rriJl 

1 J+X,J+1 



rpjl 

J ,/-1,J+1 



1 


2J + 


1 


1 




2J + 


1 


1 




2J + 


1 



jl2 T J + (j_ 

(J + l) 12 T J 



3ArpJ 



J T 



2^/j{J +1) 55 T J 
2y/j\J 



y/J(J+l)( 12 T J - 34 T J ) + 55 T J 



■J+l,J- 



_ rpjl 

1 — 1 , 7-1, .7+1 



(73) 
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The matrix elements for the five spin-dependent operators involved in Eq. ( |16[ ) in a helicity state basis, Eqs. ( J6S 
well as in \LSJM) basis, Eq. (|73|), are given in section 4 of Ref. |32|. Note that, for the amplitudes Tj^_ 1 J+1 and 
^j+i J _ 1 , we use a sign convention that differs by a factor (—1) from the one used in Ref. [ B2| . 

We consider neutron-proton scattering and determine the CM on-shell nucleon momentum p using correct relativistic 
kinematics: 

2 _ M^T lab {T lab + 2M n ) 
P ~ (M p + M n f + 2T lab M p ' ( } 



where M p — 938.2720 MeV is the proton mass, M n — 939.5653 MeV the neutron mass |26|, and Ti a b is the kinetic 
energy of the incident nucleon in the laboratory system. 

The on-shell S'-matrix is related to the on-shell T-matrix by 



with 



S J A (T lab ) = 5l>l + 2i t j J l (p, p) , (75) 

r'Aip^^-^pT^p). (76) 

For an uncoupled partial wave, the phase shifts 8j S {Ti ab ) parametrizes the partial- wave S-matrix in the form 

SfATm) = V J /(T lab ) e ^i s m^ , (77) 

implying 

t^9S JS (T y 2Rer/|(p,p) 

tan 26 j (T lab ) ——-^—^ . (78) 

The real parameter rjj s (Ti ab ), which is given by 

n J j S {T lab ) = \S J j S j(T lab )\ , (79) 

tells us to what extend unitarity is observed (ideally, it should be unity). 

For coupled partial waves, we use the parametrization introduced by Stapp et al. p3[ | (commonly known as 'bar' 
phase shifts, but we denote them simply by 5± and ej), 




\ ( cos2ej isin2ej 
(rii)ie iS + I V«sin2ej cos2ej 




(80) 



where the subscript '+' stands for 'J+1' and '— ' for 'J — V and where the superscript S = 1 as well as the argument 
Ti ab are suppressed. The explicit formulae for the resulting phase parameters are, 



Im (5^ ± /cos2 £j ) 



tan 251 = t j t^ l , -4- , (81) 

± Re (5 , ^ ± /cos2e J ) V ' 

-i S J 

tan2ej^^ + ~ , (82) 

lsi + sl_ 



vi 



S J 



±± 



cos2e . 



(83) 



The parameters S± and t]± are always real, while the mixing parameter ej is real if rj± = 1 and complex otherwise. 

We note that since the T-matrix is calculated perturbatively [cf. Eq. (p7\)] , unitarity is (slightly) violated. Through 
the parameter n^ 3 , the above formalism provides precise information on the violation of unitarity. It turns out that 
for the cases considered in this paper (namely partial waves with L > 3 and Ti ab < 300 MeV) the violation of unitarity 
is, generally, in the order of 1% or less. 

There exists an alternative method of calculating phase shifts for which unitarity is perfectly observed. In this 
method — known as the i'C-matrix approach — one identifies the real part of the amplitude V with the if-matrix. For 
an uncoupled partial-wave, the S'-matrix element, Sl, is defined in terms of the (real) i^-matrix element, kl, by 
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S L (T lab )= ] + iKL( f> P \ (84) 

1 - IK L (J), p) 



which guarantees perfect unitarity and yields the phase shift 



tan S L (T lab ) = KL (p,p) = -^^f p K L (p,p), (85) 



with Kl(j>,p) = ReVi(p,p). Combining Eqs. ([75|) and (p4|), one can write down the T-matrix element, Tl, that is 
equivalent to a given if -matrix element, kl, 

K L (p,p)+iK 2 L ( P ,p) 

TLiP ' P) = l + «i(p,p) ■ (86) 

Obviously, this T-matrix includes higher orders of K (and, thus, of V) such that consistent power counting is destroyed. 

The bottom line is that there is no perfect way of calculating phase shifts for a perturbative amplitude. Either 
one includes contributions strictly to a certain order, but violates unitarity, or one satisfies unitarity, but includes 
implicitly contributions beyond the intended order. To obtain an idea of what uncertainty this dilemma creates, we 
have calculated all phase shifts presented below both ways: using the T-matrix and A"-matrix approach. We found 
that the difference between the phase shifts due to the two different methods is smaller than 0.1 deg in F and G 
waves and smaller than 0.01 deg in H waves, at T a b = 300 MeV (and less at lower energies). Because of this small 
difference, we have confidence in our phase shift calculations. All results presented below have been obtained using 
the T-matrix approach, Eqs. (|75[) - (p3|) ■ 



C. Results 



For the T-matrix given in Eq. (67), we calculate phase shifts for partial waves with L > 3 and T/ a & < 300 MeV. 
At order four in small momenta, partial waves with L > 3 do not receive any contributions from contact interactions 
and, thus, the non-polynomial pion contributions uniquely predict the F and higher partial waves. The parameters 
used in our calculations are shown in Table I. In general, we use average masses for nucleon and pion, Mn and m^, as 
given in Table I. There are, however, two exceptions from this rule. For the evaluation of the CM on-shell momentum, 
p, we apply correct relativistic kinematics, Eq. (f74|), which involves the correct and precise values for the proton and 
neutron masses. For OPE, we use the charge-dependent expression, Eq. (|64|), which employes the correct and precise 
values for the charged and neutral pion masses. 

Many determinations of the LECs, Cj and di, can be found in the literature. The most reliable way to determine 
the LECs from empirical ttN information is to extract them from the irN amplitude inside the Mandelstam triangle 
(unphysical region) which can be constructed with the help of dispersion relations from empirical irN data. This 
method was used by Biittiker and Meifiner p7| . Unfortunately, the values for C2 and all di parameters obtained 
in Ref. |^7| carry uncertainties, so large that the values are useless. Therefore, in Table I, only c±, C3, and C4 are 
from Ref. p7| , while the other LECs are taken from Ref. |22j where the ttN amplitude in the physical region was 
considered. To establish a link between ttN and NN, we apply the values from the above determinations in our NN 
calculations. In general, we use the mean values; the only exception is C3, where we choose a value that is, in terms 
of magnitude, about one standard deviation below the one from Ref. |27). With the exception of C3, our results do 
not depend sensitively on variations of the LECs within the quoted uncertainties. 

In Figs. Uj^, we show the phase-shift predictions for neutron-proton scattering in F, G, and H waves for laboratory 
kinetic energies below 300 MeV. The orders displayed are defined as follows: 

• Leading order (LO) is just OPE, Eq. (||). 

• Next-to-leading order (NLO) is OPE plus iterated OPE, Eq. (|65|), plus the contributions of Sec. III. A (order 
two), Eqs. (0) and @. 

• Next-to-next-to-leading order (denoted by N2LO in the figures) consists of NLO plus the contributions of 
Sec. IILB (order three), Eqs. 

• Next-to-next-to-next-to-leading order (denoted by N3LO in the figures) consists of N2LO plus the contributions 
of Sec. III.C (order four), Eqs. (|9|), Q, ©-(HI), ©-©, and To this order, the phase shifts 



have never been calculated before. 
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It is clearly seen in Figs. ^-|| that the leading order 2ir exchange (NLO) is a rather small contribution, insufficient to 
explain the empirical facts. In contrast, the next order (N2LO) is very large, several times NLO. This is due to the 
irirNN contact interactions proportional to the LECs Ci that are introduced by the second order Lagrangian C n ^, 
Eq. ([)]). These contacts are supposed to simulate the contributions from intermediate A-isobars and correlated 2ir 
exchange which are known to be large (see, e. g., Ref. f36|). 

All past calculations of NN phase shifts in peripheral partial waves stopped at order N2LO (or lower). This was 
very unsatisfactory, since to this order there is no indication that the chiral expansion will ever converge. The novelty 
of the present work is the calculation of phase shifts to N3LO (the details of which are shown in Appendix A). 
Comparison with N2LO reveals that at N3LO a clearly identifiable trend towards convergence emerges (Figs. 
In G (except for 3 Gs, a problem that is discussed in Appendix A) and H waves, N3LO differs very little from N2LO, 
implying that we have reached convergence. Also 1 F$ and 3 i<4 appear fully converged. In 3 F 2 and 3 F 3 , N3LO differs 
noticeably from N2LO, but the difference is much smaller than the one between N2LO and NLO. This is what we 
perceive as a trend towards convergence. 

In Figs. we conduct a comparison between the predictions from chiral one- and two-pion exchange at N3LO 
and the corresponding predictions from conventional meson theory (curve 'Bonn'). As representative for conventional 
meson theory, we choose the Bonn meson-exchange model for the NN interaction [j36f , since it contains a compre- 
hensive and thoughtfully constructed model for 27T exchange. This 2ir model includes box and crossed box diagrams 
with NN, NA, and AA intermediate states as well as direct tttt interaction in S- and P-waves (of the tttt system) 
consistent with empirical information from irN and tttt scattering. Besides this the Bonn model also includes (repul- 
sive) w-meson exchange and irreducible diagrams of tt and p exchange (which are also repulsive). In the phase shift 
predictions displayed in Figs. [?]-||, the Bonn calculation includes only the OPE and 2tt contributions from the Bonn 
model; the short-range contributions are left out to be consistent with the chiral calculation. In all waves shown (with 
the usual exception of 3 Gs), we see, in general, good agreement between N3LO and Bonn p7| . In 3 F 2 and 3 F3 above 
150 MeV and in 3 i<4 above 250 MeV the chiral model to N3LO is more attractive than the Bonn 2-7T model. Note, 
however, that the Bonn model is relativistic and, thus, includes relativistic corrections up to infinite orders. Thus, 
one may speculate that higher orders in chiral perturbation theory (%PT) may create some repulsion, moving the 
Bonn and the chiral predictions even closer together |Q . 

The 2tt exchange contribution to the NN interaction can also be derived from empirical ttN and tttt input using 
dispersion theory, which is based upon unitarity, causality (analyticity) , and crossing symmetry. The amplitude 
NN tttt is constructed from ttN —> ttN and ttN — > ttttN data using crossing properties and analytic continuation; 
this amplitude is then 'squared' to yield the NN amplitude which is related to NN by crossing symmetry [Rj|. The 
Paris group p0| pursued this path and calculated NN phase shifts in peripheral partial waves. Naively, the dispersion- 
theoretic approach is the ideal one, since it is based exclusively on empirical information. Unfortunately, in practice, 
quite a few uncertainties enter into the approach. First, there are ambiguities in the analytic continuation and, 
second, the dispersion integrals have to be cut off at a certain momentum to ensure reasonable results. In Ref. |36|, 
a thorough comparison was conducted between the predictions by the Bonn model and the Paris approach and it 
was demonstrated that the Bonn predictions always lie comfortably within the range of uncertainty of the dispersion- 
theoretic results. Therefore, there is no need to perform a separate comparison of our chiral N3LO predictions with 
dispersion theory, since it would not add anything that we cannot conclude from Figs. fj||]. 

Finally, we like to compare the predictions with the empirical phase shifts. In G (except 3 Gs) and H waves there 
is excellent agreement between the N3LO predictions and the data. On the other hand, in F waves the predictions 
above 200 MeV are, in general, too attractive. Note, however, that this is also true for the predictions by the Bonn 
7r + 27r model. In the full Bonn model, also (repulsive) to and irp exchanges are included which bring the predictions 
to agreement with the data. The exchange of a uj meson or combined Trp exchange are 37r exchanges. Three-pion 
exchange occurs first at chiral order four. It has be investigated by Kaiser || and found to be totally negligible, at 
this order. However, 3tt exchange at order five appears to be sizable [ fl0| and may have impact on F waves. Besides 
this, there is the usual short-range phenomenology. In x?T, this short-range interaction is parametrized in terms of 
four-nucleon contact terms (since heavy mesons do not have a place in that theory) . Contact terms of order six are 
effective in P-waves. In summary, the remaining small discrepancies between the N3LO predictions and the empirical 
phase shifts may be straightened out in fifth or sixth order of xPT. 

V. CONCLUSIONS AND FURTHER DISCUSSION 

We have calculated the phase shifts for peripheral partial waves (L > 3) of neutron-proton scattering at order four 
(N 3 LO) in xPT. The two most important results from this study are: 

• At N 3 LO, the chiral expansion reveals a clearly identifiable signature of convergence. 
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• There is good agreement between the N 3 LO prediction and the corresponding one from conventional meson 
theory as represented by the Bonn Full Model J3(|. 

The conclusion from the above two facts is that the chiral expansion for the NN problem is now under control. As a 
consequence, one can state with confidence that the xPT approach to the NN interaction is a valid one. 

Besides the above fundamentally important statements, our study has also some more specific implications. A 
controversial issue that has recently drawn a lot of attention Q] is the question whether the LECs extracted from 
irN are consistent with NN. After discussing dispersion theory in the previous section, one may wonder how this 
can be an issue in the year of 2002. In the early 1970's. the Stony Brook |p2] , |39| and the Paris |43|j40]| groups showed 
independently that ttN and NN are consistent, based upon dispersion-theoretic calculations. Since dispersion theory 
is a model-independent approach, the finding is of general validity. Therfore, if 30 years later a specific theory has 
problems with the consistency of ttN and AW, then that theory can only be wrong. Fortunately, we can confirm that 
%PT for ttN and NN does yield consistent results, as we will explain now in more detail. 

The reliable way to investigate this issue is to use an approach that does not contain any parameters except for the 
LECs. This is exactly true for our calculations since we do not use any cutoffs and calculate the T matrix directly up 
to a well defined order. We then vary the LECs within their one-standard deviation range from the ttN determinations 
(cf. Table I) . We find that these variations do not create any essential changes of the predicted peripheral NN phase 
shifts shown in Figs. f|||, except for C3. Thus, the focus is on c 3 . We find that C3 = —3.4 GeV" 1 is consistent with the 
empirical phase shifts as well as the results from dispersion theory and conventional meson theory as demonstrated 
in Figs. M-^. This choice for C3 is within one standard deviation of its ttN determination and, thus, the consistency of 
ttN and NN in x?T at order four is established. 

In view of the transparent and conclusive consideration presented above, it is highly disturbing to find in the 
literature very different values for C3, allegedly based upon NN. In Ref. pj| , it its claimed that the value C3 = 
—5.08 ± 0.28 GeV -1 emerges from the world pp data below 350 MeV, whereas Ref. |4l| asserts that C3 = —1.15 
GeV -1 is implied by the JVJV phase shifts. The two values differ by more than 400% which is reason for deep concern. 

In Fig [n], we show the predictions at order four for the three values for C3 under debate. We have chosen the 
3 Fi as representative example of a peripheral partial wave since it has a rather large contribution from 2ir exchange. 
Moreover, the LEC C4 is ineffective in 3 i<4 such that differences in the choices for C4 do not distort the picture in this 
partial wave. This fact makes 3 F± special for the discussion of C3. 

Figure [fl] reveals that the chiral 2n exchange depends most sensitively on C3. It is clearly seen that the Nijmegen 
choice C3 = —5.08 GeV" 1 J2l|] leads to too much attraction, while the value C3 = —1.15 GeV -1 , advocated in Ref. [pi] , 
is far too small (in terms of magnitude) since it results in an almost vanishing 2tt exchange contribution — quite in 
contrast to the empirical NN facts, the dispersion theoretic result, and the Bonn model. 

One reason for the difference between the Nijmegen value and ours could be that their analysis is conducted at N 2 LO, 
while we go to N 3 LO. However, as demonstrated in Figs. |^-^, N 3 LO is not that different from N 2 LO and, therefore, 
not the main reason for the difference. More crucial is the fact that, in the Nijmegen analysis, the chiral 2tt exchange 
potential, represented as local r-space function, is cutoff at r = 1.4 fm (i. e., it is set to zero for r < 1.4 fm) [Q. This 
cutoff suppresses the 2ir contribution, also, in peripheral waves. If the 2tt potential is suppressed by phenomenology 
then, of course, stronger values for C3 are necessary, resulting in a highly model-dependent determination of C3. For 
example, if we multiply all non- iterative 2ir contributions by exp[— (p 2n + p' 2 ™)/A 2n ] with A w 400 MeV and n = 2, 
then with C3 = —5.08 GeV -1 we obtain a good reproduction of the peripheral partial wave phase shifts. Note that 
A ps 400 MeV is roughly equivalent to a r-space cutoff of about 0.5 fm, which is not even close to the cutoff used in 
the Nijmegen analysis. In fact, the Nijmegen r-space cutoff of r — 1.4 fm is equivalent to a momentum-space cutoff 
A w TO?,- which is bound to kill the 2tt exchange contribution (which has a momentum-space range of 2m v and larger). 
To revive it, unrealistically large parameters are necessary. 

The motivation underlying the value for C3 advocated in Ref. Q], is quite different from the Nijmegen scenario. In 
Ref. C3 was adjusted to the D waves of JVJV scattering, which are notoriously too attractive. With their choice, 
c 3 = —1.15, the D waves are, indeed, about right, whereas the F waves are drastically underpredicted. This violates 
an important rule: The higher the partial, the higher the priority. The reason for this rule is that we have more trust 
in the long-range contributions to the nuclear force than in the short-range ones. The ir + 2n contributions to the 
nuclear force rule the F and higher partial waves, not the D waves. If D waves do not come out right, then one can 
think of plenty of short-range contributions to fix it. If F and higher partial waves are wrong, there is no fix. 

In summary, a realistic choice for the important LEC C3 is -3.4 GeV -1 and one may deliberately assign an uncertainty 
of ±10% to this value. Substantially different values are unrealistic as clearly demonstrated in Figure O 
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APPENDIX A: DETAILS OF ORDER-FOUR CONTRIBUTIONS TO PERIPHERAL PARTIAL WAVE 

PHASE SHIFTS 

The order four, consists of very many contributions (cf. Sec. III.C and Figs. || and Here, we show how the 
various contributions of order four impact NN phase shifts in peripheral partial waves. For this purpose, we display 
in Fig. [Io| phase shifts for four important peripheral partial waves, namely, 1 Fs, 3 i<3, 3 Fi, and 3 G5. In each frame, 
the following individual order-four contributions are shown: 

• cf graph, first row of Fig. |, Eqs. @ and @, denoted by 'c2' in Fig. [n]. 

• Ci/Mjv contributions (denoted by 'c/M'), second row of Fig. [| Eqs. (31)-(§5j). 

• l/Ml corrections ('1/M2'), row three to six of Fig. § Eqs. 

• Two-loop contributions without the terms proportional to d, ('2-L'); Fig. [| but without the solid square; 
Eqs. (ff)-©, but with all d, = 0. 

Two- loop contributions including the terms proportional to di (denoted by 'd' in Fig. |l^); Fig. S Eqs. (|55|)-(|6^) 
with the di parameters as given in Table I. 



Starting with the result at N2LO, curve (1), the individual N3LO contributions are added up successively in the order 
given in parenthesis next to each curve. The last curve in this series, curve (6), is the full N3LO result. 

The cf graph generates large attraction in all partial waves [cf. differences between curves (1) and (2) in Fig. |§. 
This attraction is compensated by repulsion from the Ci/M/v diagrams, in most partial waves; the exception is T<3 
where Ci/M^ adds more attraction [curve (3)]. The corrections [difference between curves (3) and (4)] are 

typically small. Finally, the two-loop contributions create substantial repulsion in F% and 3 G$ which brings 1 F$ into 
good agreement with the data while causing a discrepancy for 3 G$. In 3 F 3 and 3 i<4, there are large cancelations 
between the 'pure' two-loop graphs and the di terms, making the net two-loop contribution rather small. 

A pivotal role in the above game is played by Ws, Eq. (|33), from the Ci/Mjy group. This attractive term receives a 
factor nine in 1 Fs, a factor (—3) in 3 G^, and a factor one m 3 ^ and 3 F^. Thus, this contribution is very attractive in 
F3 and repulsive in 3 G§. The latter is the reason for the overcompensation of the cf graph by the Ci/Mfj contribution 
in 3 G§ which is why the final N3LO result in this partial wave comes out too repulsive. One can expect that 1/M/v 
corrections that occur at order five or six will resolve this problem. 

Before finishing this Appendix, we like to point out that the problem with the 3 Gs is not as dramatic as it may 
appear from the phase shift plots — for two reasons. First, the 3 Gs phase shifts are about one order of magnitude 
smaller than the F and most of the other G phases. Thus, in absolute terms, the discrepancies seen in 3 Gs are small. 
In a certain sense, we are looking at 'higher order noise' under a magnifying glass. Second, the G5 partial wave 
contributes 0.06 MeV to the energy per nucleon in nuclear matter, the total of which is —16 MeV. Consequently, small 
discrepancies in the reproduction of 3 Gs by a NN interaction model will have negligible influence on the microscopic 
nuclear structure predictions obtained with that model. 



[1] S. Weinberg, Phys. Lett. B 251, 288 (1990); Nucl. Phys. B363, 3 (1991). 

[2] C. Ordonez and U. van Kolck, Phys. Lett. B 291, 459 (1992). 

[3] C. Ordonez, L. Ray, and U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994); Phys. Rev. C 53, 2086 (1996). 

[4] U. van Kolck, Prog. Part. Nucl. Phys. 43, 337 (1999). 

[5] L. S. Celenza, A. Pantziris, and C. M. Shakin, Phys. Rev. C 46, 2213 (1992). 



14 



[6] C. A. da Rocha and M. R. Robilotta, Phys. Rev. C 49, 1818 (1994); ibid. 52, 531 (1995); M. R. Robilotta, Nucl. Phys. 
A595, 171 (1995); M. R. Robilotta, and C. A. da Rocha, Nucl. Phys. A615, 391 (1997); J.-L. Ballot, C. A. da Rocha, 
and M. R. Robilotta, Phys. Rev. C 57, 1574 (1998). 
[7] N. Kaiser, R. Brockmann, and W. Weise, Nucl. Phys. A625, 758 (1997). 
[8] N. Kaiser, S. Gerstendorfer, and W. Weise, Nucl. Phys. A637, 395 (1998). 
[9] N. Kaiser, Phys. Rev. C 61, 014003 (1999); ibid. 62, 024001 (2000). 
[10] N. Kaiser, Phys. Rev. C 63, 044010 (2001). 
[11] N. Kaiser, Phys. Rev. C 64, 057001 (2001). 
[12] N. Kaiser, Phys. Rev. C 65, 017001 (2002). 

[13] E. Epelbaum, W. Glockle, and U.-G. Meifiner, Nucl. Phys. A637, 107 (1998); ibid. A671, 295 (2000). 
[14] D. R. Entem and R. Machleidt, Phys. Lett. B524, 93 (2002). 

[15] D. B. Kaplan, M.J. Savage, and M.B. Wise, Nucl. Phys. B478, 629 (1996); D. B. Kaplan, Nucl. Phys. B494, 471 (1997); 

D. B. Kaplan, M.J. Savage, and M.B. Wise, Nucl. Phys. B534, 329 (1998); Phys. Lett. B424, 390 (1998). 
[16] R. J. Furnstahl, B. D. Serot, and H.-B. Tang, Nucl. Phys. A615, 441 (1997); J. V. Steel and R. J. Furnstahl, ibid. A637, 

46 (1998); R. J. Furnstahl, J. V. Steel, and N. Tirfessa, ibid. A671, 396 (2000); R. J. Furnstahl, H. W. Hammer, and N. 

Tirfessa, ibid. A689, 846 (2001). 
[17] T.-S. Park, K. Kubodera, D. P. Min, and M. Rho, Phys. Rev. C 58, 637 (1998); Nucl. Phys. A646, 83 (1999). 
[18] T. D. Cohen, Phys. Rev. C 55, 67 (1997); D. R. Phillips and T. D. Cohen, Phys. Lett. B390, 7 (1997); K. A. Scaldeferri, 

D. R. Phillips, C. W. Kao, and T. D. Cohen, Phys. Rev. C 56, 679 (1997); S. R. Beane, T. D. Cohen, and D. R. Phillips, 

Nucl. Phys. A632, 445 (1998). 

[19] G. Rupak and N. Shoresh, |nucl-th/9906077| ; P. F. Bedaque, H. W. Hammer, and U. van Kolck, Nucl. Phys. A676, 357 

(2000); S. Fleming, Th. Mehen, and I. W. Stewart, Nucl. Ph ys. A677, 313 (2 000); Phys. Rev. C 61, 044005 (2000). 
[20] S. R. Bean, P. F. Bedaque, M. J. Savage, and U. van Kolck, (tiucl-th/0104030 . 

[21] M. C. M. Rentmeester, R. G. E. Timmermans, J. L. Friar, and J. J. de Swart, Phys. Rev. Lett. 82, 4992 (1999). 
[22] N. Fettes, U.-G. Meifiner, S. Steiniger, Nucl. Phys. A640, 199 (1998). 

[23] N. Fettes, U.-G. Meifiner, M. Mojzis, and S. Steininger, Ann. Phys. (N.Y.) 283, 273 (2000); ibid. 288, 249 (2001). 

[24] V. Stoks, R. Timmermans, and J. J. de Swart, Phys. Rev. C 47, 512 (1993). 

[25] R. A. Arndt, R. L. Workman, and M. M. Pavan, Phys. Rev. C 49, 2729 (1994). 

[26] Review of Particle Physics, Eur. Phys. J. C 15, 1 (2000). 

[27] P. Biittiker and U.-G. Meifiner, Nucl. Phys. A668, 97 (2000). 

[28] V. Bernard, N. Kaiser, and U.-G. Meifiner, Int. J. Mod. Phys. E 4, 193 (1995). 

[29] G. Q. Li and R. Machleidt, Phys. Rev. C 58, 3153 (1998). 

[30] M. Jacob and G. C. Wick, Ann. Phys. (N.Y.) 7, 404 (1959). 

[31] J. Goto and S. Machida, Prog. Theor. Phys. 25, 64 (1961). 

[32] K. Erkelenz, R. Alzetta, and K. Holinde, Nucl. Phys. A176, 413 (1971); note that there is an error in equation (4.22) of 

this paper where it should read: ~W£s = 2g< ? '^ T [ J 4^~ 2 ' (0) - A J L f] and +W£ S = 2 OT '^[A^+ 2 ' (0) - A J L f } ]. 
[33] H. P. Stapp, T. J. Ypsilantis, and N. Metropolis, Phys. Rev. 105, 302 (1957). 

[34] V. G. J. Stoks, R. A. M. Klomp, M. C. M. Rentmeester, and J. J. de Swart, Phys. Rev. C 48, 792 (1993). 

[35] R. A. Arndt, I. I. Strakovsky, and R. L. Workman, SAID, Scattering Analysis Interactive Dial-in computer facility, Virginia 

Polytechnic Institute and George Washington University, solution SM99 (Summer 1999); for more information see, e. g., 

R. A. Arndt, I. I. Strakovsky, and R. L. Workman, Phys. Rev. C 50, 2731 (1994). 
[36] R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 

[37] Note that the Bonn model uses the tyNN coupling constant fl 2 ]vjv/47r = 14.4, while for chiral pion exchanges we apply 

Ann 1^ = 13-6 [cf. footnote a ) of Table I]. This is the main reason for the light "discrepancies" between N3LO and Bonn 

that seem to show up in some of the H waves (Fig. |^). 
[38] In fact, preliminary calculations, which take an important class of diagrams of order five into account, indicate that the 

N 4 LO contribution may prevailingly be repulsive (N. Kaiser, private communication). 
[39] G. E. Brown and A. D. Jackson, The Nucleon-Nucleon Interaction, (North-Holland, Amsterdam, 1976). 
[40] R. Vinh Mau, in: Mesons in Nuclei, Vol. I, M. Rho and D. H. Wilkinson, eds. (North-Holland, Amsterdam, 1979) p. 151; 

M. Lacombe, B. Loiseau, J. M. Richard, R. Vinh Mau, J. Cote, P. Pires, and R. de Tourreil, Phys. Rev. C 21, 861 (1980). 
[41] E. Epelbaum, A. Nogga, W. Glockle, H. K amada, U.-G. Meifi ner, and H. Witala, Few-Nucleon Systems with Two-Nucleon 

Forces from Chiral Effective Field Theory, rmcl-th/0201064. 
[42] M. Chemtob, J. W. Durso, and D. O. Riska, Nucl. Phys. B38, 141 (1972); A. D. Jackson, D. O. Riska, and B. Verwest, 

Nucl. Phys. A249, 397 (1975). 

[43] R. Vinh Mau, J. M. Richard, B. Loiseau, M. Lacombe, and W. M. Cottingham, Phys. Lett. B44, 1 (1973); M. Lacombe, 
B. Loiseau, J. M. Richard, R. Vinh Mau, P. Pires, and R. de Tourreil, Phys. Rev. D 12, 1495 (1975). 

[44] We note that in the Nijmegen analysis [ElJ C4 = 4.7 ± 0.7 Gev -1 emerges which also differs from what we use and what is 
obtained in nN analysis (cf. Table I). However, this difference in a has very little impact on NN peripheral partial waves 
and does not make possible stronger values for C3. 



15 



TABLE I. 


Parameters used in our calculations. The LECs Ci and c 


i are in units of GeV 1 and GeV ^, respectively. 




Our choice 


Empirical 


M N 


938.9182 MeV 




m„ 


138.039 MeV 




9A 


1.29 


1.29 ±0.01° 


fn 


92.4 MeV 


92.4 ± 0.3 MeV 6 


Cl 


-0.81 


-0.81 ± 0.15 c 


C2 


3.28 


3.28 ± 0.23 d 


C3 


-3.40 


-4.69±1.34 c 


C 4 


3.40 


3.40 ± 0.04 c 


di + d 2 


3.06 


3.06±0.21 d 


d 3 


-3.27 


-3.27 ± 0.73 d 


d 5 _ 


0.45 


0.45 ± 0.42 d 


dl4 — dl5 


-5.65 


-5.65±0.41 d 



"Using g^ NN /^ — 13.63 ± 0.20 |3J25fl and applying the Goldberger-Treiman relation, qa = g-nNN fn/M^. 



b ReL [[26|. 
c Table 1, Fit 1 of Ref. 
d Table 2, Fit 1 of Ref. 
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FIG. 1. Two-pion exchange contributions to the NN interaction at order two and three in small momenta. Solid lines 
represent nucleons and dashed lines pions. Small dots denote vertices from the leading order ttN Lagrangian £^r, Eq. 

Eq. (|l|). Symbols with 
Eqs. (0). Only a few 



(2) 



Large solid dots are vertices proportional to the LECs c,: from the second order Lagrangian 

an open circles are relativistic 1/Mjv corrections contained in the second order Lagrangian 
representative examples of 1/Mjv corrections are shown and not all. 
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16 
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(N 3 LO) 



4 , 



(®) 1 <ii (i) <i 




FIG. 2. One-loop 27r-exchange contributions to the NN interaction at order four. Basic notation as in Fig. |l|. Symbols with 
a large solid dot and an open circle denote 1/Mjv corrections of vertices proportional to C;. Symbols with two open circles mark 
relativistic 1/Mjf corrections. Both corrections are part of the third order Lagrangian Eq. (|l4|). Representative examples 

for all types of one-loop graphs that occur at this order are shown. 




FIG. 3. Two-loop 27r-exchange contributions at order four. Basic notation as in Fig. [y. The grey disc stands for all one-loop 
ttN graphs some of which are shown in the lower part of the figure. The solid square represents vertices proportional to the 
LECs di which are introduced by the third order Lagrangian C^ N , Eq. (jis]). More explanations are given in the text. 
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FIG. 4. F-wave phase shifts of neutron-proton scattering for laboratory kinetic energies below 300 MeV. We show the predic- 
tions from chiral pion exchange to leading order (LO), next-to-leading order (NLO), next-to-next-to- leading order (N2LO), and 
next-to-next-to-next-to-leading order (N3LO). The solid dots and open circles are the results from the Nijmegen multi-energy 
np phase shift analysis pm and the VPI single-energy np analysis SM99 [B5|, respectively. 
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FIG. 5. Same as Fig. W, but for G-waves. 
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FIG. 6. Same as Fig. ^, but for Jf-waves. 
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FIG. 7. _F-wave phase shifts of neutron-proton scattering for laboratory kinetic energies below 300 MeV. We show the results 
from one-pion-exchange (OPE), and one- plus two-pion exchange as predicted by xPT at next-to-next-to-next-to- leading order 
(N3LO) and by the Bonn Full Model |Q (Bonn). Empirical phase shifts (solid dots and open circles) as in Fig. ^. 
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FIG. 8. Same as Fig. [j], but for G-waves. 
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FIG. 9. Same as Fig. ^, but for _ff-waves. 



20 



-1 



\ 



■a 



<n 
re 



\ 



+C/M 
/M2 



. (3) +c 
sS- : - ^ (5) +2-L 

o 



(1) N2L0 
X (6) +d 



O 

_i_ 



100 200 300 
Lab. Energy (MeV) 



400 



0) 

5- 



in 

0) 

« -1 




. (4) +1/M2 

5)+2-L 
3 +C/M 
1 N2L0 



5 

• s 



100 200 300 400 
Lab. Energy (MeV) 



C/3 

0) 

CO 
re 
.n 
Q. 




100 200 300 400 
Lab. Energy (MeV) 



0.25 



3Q /(2 ) + c2 



/ / 
/ 

y 



(1) N2LO 




v*-Q - (3) +C/M 

^0^0..0(4) +1/M2 

N 

\ 

^ (5) +2-L 



100 200 300 400 
Lab. Energy (MeV) 



FIG. 10. The effect of individual order-four contributions on the neutron-proton phase shifts in some selected peripheral 
partial waves. The individual contributions are added up successively in the order given in parenthesis next to each curve. 
Curve (1) is N2LO and curve (6) is the complete N3LO. For further explanations, see Appendix A. Empirical phase shifts (solid 
dots and open circles) as in Fig. ^. 
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FIG. 11. One and two-pion exchange contributions at order four to the 3 f4 phase shifts for three different choices of the LEC 
C3. The numbers given next to the curves denote the values for C3 in units of GeV~ used for the respective curves (all other 
parameters as in Table I). For comparison, we also show the OPE contributrion (OPE) and the result from n + 2n exchange 
of the Bonn model (Bonn). Empirical phase shifts (solid dots and open circles) as in Fig. m 
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